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CMB Cosmological Parameter Estimation: Methods and Current Results 
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The majority of present efforts to constrain cosmological parameters with cosmic microwave 
background (CMB) anisotropy data employ approximate likelihood functions, the time con- 
suming nature of a complete analysis being a major obstacle. We have performed a full 
(unapproximated) likelihood analysis on several experiments that allows us to examine the 
various assumptions made in these approximate methods and to evaluate their performance. 
Our results indicate that care should be taken when using such approaches. With an improved 
approximate method, we present some constraints on cosmological parameters using the entire 
present CMB data set. 



1 Introduction 

Since their first detection by COBE, the CMB temperature fluctuations have become an es- 
sential tool for constraining cosmological parameters. The number of experiments measuring 
anisotropics at different scales is now close to 20. With time, the observations become more 
accurate and their angular resolution improves, increasing the number of pixels. This evolution 
towards higher precision measurements of CMB fluctuations inhibits complete (unapproximated) 
likelihood analyses over large sets of models (or free parameters)[|. Approximate methods have 
thus been proposed in order to explore as much parameter space as possible with the complete 
CMB data set. Easy to use and fast, minimisation has been widely applied. The approximate 
nature of this method resides in two assumptions: (i) the one dimensional bandpower distribu- 
tion is Gaussian; (ii) all pertinant information in a map is contained in the flat-band power. 
The latter is also assumed by other approximations, such as proposed by Bond et al. (1998), 
Wandelt et al. (1998) and Bartlett et al. (1999). Focussing on a small number of experiments 
(COBE, SASKATOON and MAX), we have performed a complete likelihood analysis over a 
reasonably representative set of (Inflationary) models. This allows us to compare our results 
with those from a minimisation and to quantitatively examine these two assumptions. 



^note that the computation time grows like A''pi, 
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Figure 1: True likelihood (solid line), 2-tailed Gaussian (dot-dot-dashed) and our approximation (Bartlett et al. 

1999, dot-dashed) for the MAX ID field. 

2 Testing the assumptions 

In an inflationary scenario, the sky pixels are Gaussian distributed, and so also the a'l^s in a 
spherical harmonic decomposition. But what is used in minimisation is the in-band power 
of the fluctuations; this corresponds to the square root of the temperature fluctuation variance 
and is not a Gaussian distributed quantity. This motivated us to look more carefully at the first 
assumption. Figure 1 shows an example where a 2-tailed Gaussian (with asymmetric errors; 
dot-dot-dashed line) is a rather bad approximation to the true likelihood curve (solid line). 
Even if the Gaussian appears relatively faithful near the maximum, it rapidly deviates from 
the likelihood function with the distance from the peak: the is thus very (overly) sensitive 
to the presence of outliers (which will be more common than it expects). By comparison, 
the approximation developed in Bartlett et al. (1999) reproduces quite well the likelihood 
function[[ More accurate (relative to a full likelihood treatment) parameter estimation can thus 
be obtained by using such approximate likelihood forms, or simply by direct interpolation of the 
one-dimensional in-band likelihood function, when available. 

The second assumption supposes that the flat -band power estimate contains all the infor- 
mation in a map. To check this assumption, we computed the flat-band power and its error 
bars for several bins of MAX and Saskatoon and will now compare the constraints suggested 
by these estimates to those from a full likelihood analysis. This is illustrated in Figure 2. By 
definition, models passing within the error bars are considered consistent with the flat-band 
estimates at the "one sigma confidence level". This may be compared to the confidence level 
assigned by the full likelihood analysis: focussing on just the MAX PH point (at the extreme left 
in Figure 2-left), the likelihood analysis excludes the model plotted as the dotted line at more 
than 68%, but accepts the model shown as the dashed line, all in flagrent contradiction with the 
flat-band estimate. This can be understood by looking at Figure 2-right. Here, we show the 
results of a likelihood analysis over a family of spectra modeled by their in-band power ((5T/v) 
and slope (m): 5T{£) = 5TN.{£/£eff)^- Figure 2~right displays likelihood contours in the plane 



''Many examples of the comparison between a Ga ussian, the likelihood function and our approximation may 
be seen at http://webast.ast.obs-mip.fr/cosmo/CMB 




{5T]\f,m) for MAX PH. We clearly see that MAX PH prefers spectra with positive slope and 
lower in-band power than the flat -band estimate (corresponding to the vertical line at m=0). 
The latter considers only the information available along the vertical line in Figure 2-right and 
looses signal contained in the rest of the plane. The symbols position the two models plotted in 
Figure 2-left, explaining the disagreement between the two analysis methods. The conclusion 
is that, at least in this example, the information contained in the original pixel set cannot be 
reduced to a flat-band estimate. In the case of MAX PH, any approximation using only the 
flat-band power will never result in the same constraints as a full likelihood analysis. 

Incorporating all the information available in a map (in-band power, slope, ...)[] will give 
more realistic constraints than possible with just flat-band estimates. A simple and efficient 
way to impliment this is by interpolation the likelihood surface over JT/v and m for each bin 
(instead of the one-dimensional curve for the flat-band power) . 

3 Current Results 

We took into account the above remarks in an attempt to obtain a more complete parameter 
estimation using the entire CMB data set and a large set of Inflationary models. In practice, 
we used our approximation for those experiments for which only the flat-band power estimate 
was given; interpolated the flat-band likelihood curve when it was available; and interpolated 
the likelihood surfaces (power, slope) of Saskatoon and MAX. Our analysis corresponds to 
70 different data points representing 20 experiments, and more than 10 million models. The 
compute time required was 350 cpu-hours on a DEC workstation - equivalent to a simple 
minimisatiorl[|. 

Figures 3 show our principle results in the framework of open Inflationary models with a 
cosmological constant. The constraints are given as 2-dimensional (approximated) likelihood 
contours at roughly 68 and 95% confidence. In each case, the other investigated parameters 

'^Depending on the signal-to-noise, one could imagine the need for additional parameters. 
"^For comparison, the estimated time for a full likelihood treatment with the same data set and models would 
be ~ 10^ cpu-hours 




Figure 3; LEFT: Constraints in {fltot,Xo) plane obtained with our method. RIGHT: Constraints in in the 

{Qtot,n) plane. 



{Hq, Q, n, r/io) are marginalised by projection. A large number of the models are excluded in the 
plane (Jlfc, Aq)- Prefered models have very low curvature, although no the cosmological constant 
is relatively unconstrained. Figure 3-right is a nice summary of the status of the simplest of 
Inflationary predictions, namely that n ~ 1 and 0,k = 0. These values are indeed favored by 
current CMB observations. 



4 Conclusion 



We have seen how in some cases a Gaussian is a bad representation of the one-dimensional 
in-band power likelihood function. To improve on methods, we therefore recommend the use 
of other, more appropriate approximations, as proposed in the literature, or direct interpolation 
of the exact likelihood curve, when available. Even this may not be entirely satisfactory be- 
cause flat -band power estimates do not always sufficiently represent a CMB data set. We have 
demonstrated that some of MAX and Saskatoon data (bins) prefer non-zero spectral slopes and 
different in-band powers. In these cases, the incorporation of additional information, for exam- 
ple, the in-band spectral slope, should lead to a better reproduction of the actual constraints. 
Taking these remarques into account, have also seen how the first generation experiments favor 
the simple Inflationary scenario predictions: curvature close to zero and spectral index close to 
one. 
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